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Abstract 

We present simulations performed with REAS2, a new Monte Carlo code for the cal- 
culation of geosynchrotron radio emission from extensive air showers. The code uses 
thoroughly tested time-domain radio emission routines in conjunction with a real- 
istic air shower model based on per-shower multi-dimensional CORS IK A- generated 
histograms. We assess in detail how the transition from simpler, parametrised, to 
realistic, CORSIKA-based particle distributions affects the predicted radio emission 
from a typical 10^^ eV air shower. The effects of eliminating a previously needed free 
parameter and adopting realistic electron to positron ratios are also discussed. Com- 
pared with earlier calculations based on parametrised showers, REAS2 simulations 
predict slightly weaker and in some cases narrower pulses. In addition, a pronounced 
east- west versus north-south asymmetry arises in the emission pattern, and the ra- 
dio pulses become generally unipolar. Finally, we demonstrate how REAS2 can be 
used to study radio pulse shapes and their relation to shower characteristics such 
as the longitudinal air shower development. 
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1 Introduction 



In the last few years, radio detection of cosmic ray air showers has once again 
become a very active field of research. A number of new experiments such 
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as LOPES [H; H 3 3 and CODALEMA 0; Q have been established. These 
experimental activities are accompanied by theoretical studies of the emission 
processes to build the necessary basis for the interpretation of experimental 
data. 

There have already been a number of pioneering studies of radio emission from 
extensive air showers (EAS) in the context of geomagnetic and Cerenkov-like 
mechanisms (see for a review). A new perspective on radio emission from 
EAS was introduced with the proposal to interpret the dominant geomagnetic 
contribution as synchrotron radiation from secondary shower electrons and 
positrons gyrating in the earth's magnetic field j^, hereafter called "geosyn- 
chrotron radiation" . 

Analytical calculations along that line have been presented in and, in as- 



sociation with the LOPES project, in [10|. In the model of Huege & Falcke 



[10| , the longitudinal evolution of the air shower as well as the lateral, arrival 
time and energy distributions of the particles are described with well-known 
analytical parametrisations. The radio emission is calculated in the frequency- 
domain. Analytical calculations like these help in understanding complex co- 
herence effects arising from air shower structures on length scales comparable 
to the observing wavelength. On the other hand, analytical calculation in the 
frequency-domain necessitate certain approximations such as adopting infinite 
particle tracks and the adoption of the Fraunhofer limit. A recent study [ll| 
takes into account the effect of semi-infinite particle tracks. A general draw- 
back of analytical calculations, however, remains: they are very limited in the 
description of the extremely complex air shower characteristics. 

A natural choice to tackle this complexity is to revert to Monte Carlo simu- 



lations. In [12] we have described a Monte Carlo code for the calculation of 
geosynchrotron radio emission from air showers (REASl). The code calcu- 
lates the emission in the time-domain. This approach has several advantages: 
no far-field approximations have to be made, the finiteness of the particle 
trajectories is taken into account in a natural way, and polarisation charac- 
teristics of the emission can be easily calculated. In order to allow a good 



comparison with the results of [10|, the air shower model was kept identical. 



i.e., the air shower evolution and particle distributions were described with 



the same analytical parametrisations as in [10|. This code allowed extensive 



analyses of the relationship between air shower parameters and the associated 



radio emission, the results of which have been published in [13 



Other efforts to simulate radio emission from cosmic ray air showers with 
Monte Carlo techniques have recently been undertaken. Essentially the same 



formalism that was used in 12(] has been incorporated into the AIRES 14] 
code for extensive air shower simulations. Some first resu lts p resented in 15] 



showed qualitatively similar results to those published in [13|] . In a somewhat 
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different approach, the EGSnrc [16| code has been adapted to make frequency- 
domain radio calculations taking into account the refractive index profile of the 
atmosphere, i.e., incorporating both the geomagnetic as well as an Askaryan- 
type Cerenkov radiation mechanism. Initial results at energies around 10 TeV 
have been presented in [l7|. 

This article describes the next step (REAS2) in the development of our geosyn- 
chrotron radiation model. The well-tested radio emission part of the Monte 
Carlo code is essentially kept the same as in REASl, but the analytical descrip- 
tion of the air shower characteristics is replaced by detailed multi-dimensional 



histograms obtained from state-of-the-art CORSIKA [18j simulations. The 
much more realistic air shower model allows calculations of the radio emission 
with much higher precision. At the same time, the code allows us to make 
a very gradual transition from the earlier, well-understood and thoroughly 
tested simulations, providing insights in the importance of various aspects of 
particle distributions in EAS for the radio signal. 

In this article, we first discuss the strategy of Monte Carlo simulations with 
REAS2 and then analyse in detail how the improvement of the underlying air 
shower model, the elimination of the previously needed track length parameter, 
and the adoption of the true electron to positron ratio change the predicted 
radio emission. In the second part we illustrate how the new simulations can 
be used to analyse different regimes identifiable in the particle distributions 
and their contribution to the radio pulses, before concluding with an outlook 
on the future development. 



2 CORSIKA-based Monte Carlo simulations 



With REAS2, the process of simulating radio emission from extensive air show- 
ers is separated into two steps. First, the air shower is simulated as usual 
with CORSIKA, adopting the desired interaction models and simulation pa- 
rameters. A tailor-made interface code is used to collect the relevant particle 
information in histograms during the air shower simulation. At the end of the 
simulation, the histograms are saved to disk in a compact data file based on 
ROOT data structures. 

The following particle information is collected for electrons and positrons sep- 
arately in (usually) 50 slices equidistantly distributed in slant atmospheric 
depth over the air shower evolution (from first interaction to observer level): 

• one three-dimensional histogram of 

(1) particle arrival time relative to that of an imaginary primary particle 
propagating with the speed of light from the point of first interaction 
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(2) lateral distance of the particle from the shower core 

(3) particle energy 

• and one three-dimensional histogram of 

(1) angle of the particle momentum to the shower axis 

(2) angle of the particle momentum to the (radial) outward direction 

(3) particle energy 

With the first set of histograms we have access to a four-dimensional histogram 
of atmospheric depth, particle arrival time, lateral distance and energy. In 
addition, we have information on the particle momentum angles to the shower 
axis and outward direction as a function of atmospheric depth and particle 
energy. For atmospheric depth values between the 50 slices, the histograms are 
interpolated appropriately. To improve the sampling of the shower evolution, 
the shower size longitudinal profile is taken from CORSIKA directly (sampled 
usually every 5 g cm~^). 



Naturally, one complete higher- dimensional histogram would be superior to 
two separate lower-dimensional histograms. The approximation made by sep- 
arating the information into the two above-mentioned histograms is, however, 
motivated by air shower physics and allows us to describe the air shower 
properties well with only moderate memory usage. Regarding the angular dis- 
tribution of particle momenta, we chose to retain its correlation with particle 
energy rather than that with lateral distance. The correlation with particle 
energy is very strong as illustrated in [l9|. Fig. 11. An additional correlation 
with lateral distance certainly exists, but is neglected here. The radio signal 
is dominated by the bulk of the particles close to the shower axis, whereas 
particles at large lateral distances contribute only weakly (cf. Sec. 7.2). For 
these particles close to the shower axis, however, the correlation with lateral 
distance can only be weak; otherwise, the angular distributions plotted in Fig. 
11 of 191] would show a stronger evolution with shower age. The current ap- 
proximation gives a good representation of rotationally symmetric showers. 
Possible asymmetry effects due to the geomagnetic field (such as a systematic 
shift of electrons with respect to positrons) are not considered, but will be 
included at a later time. In any case, the histograms retain the systematic 
outward drift of particles from the shower core. 



At energies of about 10^^ eV and higher, the thinning option in CORSIKA is 
used to keep computation times acceptable. As we are only interested in the 
histogrammed particle distributions, the requirements on the quality of the 
thinned simulations is only moderate. Optimum 10~^ thinning 20] yields high- 
quality simulations clearly detailed enough for the calculation of histogrammed 
particle distributions and was therefore used in all simulations presented here. 



In a second step, these histograms are imported in the REAS2 code, where 
they are used to recreate particles which follow the given distributions and 
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then calculate the radio emission from the deflection of the particles in the 
earth's magnetic fleld. The well-understood and thoroughly tested radio emis- 
sion calculation is essentially unchanged from our earlier code (for details see 

H). 

This two-step approach has a number of advantages. In particular, it allows 
us to make a very gradual transition from fully parametrised air showers to 
fully histogrammed air showers by switching on the histogrammed distribu- 
tions individually one after the other and evaluating the changes arising in the 
radio signal. The fact that one does not have to re-run the complete air shower 
simulation each time one wants to calculate the radio emission in a different 
conflguration (e.g., with parametrised vs. histogrammed distributions) is also 
very helpful in keeping the computation time low. A reduction of computa- 
tion times has also been achieved by the implementation of a new time-grid 
algorithm and an algorithm to automatically optimise the time window and 
time resolution for a speciflc radio pulse calculation as a function of lateral 
observer distance from the shower centre. On a current standard PC (2 GHz 
class, 512 MB memory), the time distribution radio pulse calculation (without 
the CORSIKA run) takes approximately 36-48 hours for 200 observer positions 
(antennas). 



3 Effects of the improved distributions 



We investigate the effects arising due to the transition from parametrised dis- 
tributions to CORSlKA-based histograms for the example of a typical 10^^ eV 
vertical proton-induced air shower with shower maximum at ~ 640 g cm~^. A 
magnetic fleld of 0.48 Gauss with an inclination angle of 64.7°, as appropriate 
for the LOPES-location in Karlsruhe, Germany, is adopted for all simulations. 



All CORSIKA simulations were carried out with the QGSJETII.03 |2l|; l22 



high-energy and UrQMDl.3.1 [23|;|2J| low-energy interaction models. Electro 



magnetic cascades were treated with EGS4 [25[ down to a threshold of 400 keV 
kinetic energy, where the particles are no longer relativistic and therefore be- 
come unimportant for the geosjTichrotron mechanism. 

To illustrate the changes arising in the transition from the parametrised to 
the histogrammed air shower model, we investigate the particle distributions 
individually, switching from parametrised to histogrammed distributions one 
after the other in a series of flve steps. In the diagrams, we mark distributions 
which are parametrised with a "0" and distributions which are histogrammed 
with a "1" inside a 5-digit key. The order of the digits in the key corresponds 
to the order in which the histogrammed distributions replace the parametrised 
ones in the subsections 3.1 to 3.5. 
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In each of the five steps we look at the unhmited-bandwidth radio pulses at a 
close distance of 75 m and at a larger distance of 525 m to the north from the 
shower centre. The pulses at 525 m distance are amplified by a factor of 20 
for a clearer illustration. These two distances represent different regimes for 
experimental measurements. Closer distances are particularly important for 
LOPES, while larger distances are relevant for radio detection on large scales 
as currently planned for the Pierre Auger Observatory. 

In a sixth step, we then investigate how the transition to the histogrammed 
particle distributions changes the pulses seen by observers 75 m and 525 m to 
the west from the shower centre. 



3.1 Air shower longitudinal evolution 



The air shower lon gitu dinal development in REASl was parametrised with 
a Greisen function [lO|; |26| which was originally designed to describe purely 
electromagnetic showers. The position of the shower maximum in REASl was 
a user-defined parameter. Using this parameter one could mimic the effect of 
different primary particle masses, yet only in a crude way. Shower-to-shower 
fluctuations were not taken into account. 



Air showers simulated with CORSIKA naturally take into account differences 
in the shower evolution related to the primary particle mass. Also, shower-to- 
shower fluctuations and the influence of different high- and low-energy interac- 
tion models are automatically accounted for. For the analyses presented in this 
article, we simulated a number of vertical 10^^ eV proton-induced CORSIKA 
showers and then selected one with a typical proflle and a shower maximum 
close to the expected average depth of ~ 640 g cm~^. 
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Fig. 1. Particle (sum of e"^ and e~) in- 
jection profiles in parametrised and his- 
togrammed versions. 
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Fig. 2. Radio pulses for parametrised 
and histogrammed air shower evolu- 
tions. 



In flgure 1 we compare the parametrised and CORSIKA-based particle injec- 
tion proflles, which are closely related to the corresponding air shower lon- 
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gitudinal profiles via the track-length parameter A (cf. Sec. 4; throughout 
section 3, A is fixed to 36.7 g cm~^ with an exponential distribution of par- 
ticle track lengths). The parametrised profile and the profile of the selected 
CORSIKA shower have virtually identical integrated particle numbers, but 
the CORSIKA-based profile is slightly wider with a lower particle number at 
shower maximum. 

Figure 2 shows the radio pulse^l] corresponding to a fully parametrised shower 
(00000) and a shower with histogrammed air shower evolution profile, but all 
other distributions parametrised (10000). The changes to the pulse shapes are 
very small. 



3.2 Particle energies 

We now additionally switch the particle energy distribution from parametrised 
to histogrammed. Please note that the histogrammed energy distribution changes 
as a function of atmospheric depth, while the parametrised distribution was 
used globally throughout the shower. 

In figure 3 we compare the parametrised particle energy distribution with the 
histogrammed particle energy distributions as obtained from CORSIKA. The 
histogrammed energy distributions deviate significantly from the parametri- 
sation. In particular, there are more particles at lower energies and there is no 
peak in the distribution at medium Lorentz factors. 

Figure 4 shows the effect of switching on the histogrammed energy distribu- 
tion. The changes are only slight. The relative up-weighting of low energies, 
and, correspondingly, down-weighting of higher energies is mainly visible in a 
reduced pulse height in the shower centre, where high-energy particles with 
narrow beaming cones contribute strongly (cf. Sec. 7). The disappearance of 
the prominent peak in the energy distribution around Lorentz factors of 60 is 
important for the pulse shapes along the east-west axis from the shower centre 
(cf. Sec. 3.6). 



Throughout this article we show radio pulses corresponding to the east-west po- 
larisation component of the emission. This allows us to retain the sign of the electric 
field while its strength is virtually identical to the total field strength (in case of 
vertical showers). Positive field strengths denote electric fields pointing to the west. 
The zero-point in time is arbitrarily set to the start of the 75 m pulses. 
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Fig. 3. Distribution of particle energies 
in parametrised and histogrammed ver- 
sions (averaged over all other distribu- 
tions, normalised to a total of 10^ par- 
ticles) . 
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Fig. 4. Radio pulses for fully 

parametrised showers compared 

with those for showers with his- 
togrammed shower evolution and 
energy distribution. 



3.3 Particle lateral distribution 



In this step, we additionally switch on the histogrammed lateral particle 
distribution. The parametrised distribution already accounted for changes 
with atmospheric depth through the variation of the Moliere radius. The his- 
togrammed distribution takes into account both the development with atmo- 
spheric depth and the correlation with the particle energy distribution. 



In figure 5 we compare the lateral distributions of shower particles in the 
parametrised and histogrammed versionsH] The parametrisation used is an 



NKG-function as given in 26|; |27l|. The histogrammed distributions have a 
higher particle density in the medium- distance region, but a lower density in 
the extreme centre. 



Figure 6 shows the corresponding radio pulses. The pulse heights in the shower 
centre are slightly amplified with respect to Fig. 4, because the histogrammed 
lateral distribution is characterised by more high-energy particles near the 
shower core. 



^ The "zig-zag" structure present in the histogrammed distributions shown 
throughout this article is an artifact of the adopted logarithmic histogram binning. 
As the logarithmic bin boundaries have no regular structure on a linear scale, this 
is unproblematic for the calculation of the radio signal. In most cases logarithmic 
binning gives a better representation of the particle distributions than using linear 
intervals. 
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Fig. 5. Lateral particle distribution in 
parametrised and histogrammed ver- 
sions (averaged over all other distribu- 
tions, normalised to a total of 10^ par- 
ticles) . 
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Fig. 6. Radio pulses for fully 
parametrised showers compared with 
those for showers with histogrammed 
shower evolution, energy and lateral 
distributions. 



3.4 Particle arrival times 



We now additionally switch the particle arrival time distribution from parametrised 
to histogrammed. The parametrised distribution was only correlated with ra- 
dial distance from the shower axis. In particular, it was assumed to be the same 
over the full shower evolution. In contrast, the histogrammed arrival time dis- 
tribution is now fully correlated with atmospheric depth, lateral distance and 
particle energy. 

The arrival time distribution has strong influence on the shape (and thus 
frequency spectrum) of the radio pulses (for a discussion, see [lO[). Figure 7 
shows the parametrisation of the particle arrival times used so far in compar- 
ison with the distributions obtained from CORSIKA. The CORSIKA-based 
overall arrival time distributions are much narrower. This was to be expected 



because the parametrised distribution [10| was derived from measurements 



28l | that are very difficult if not impossible to perform close to the shower 
axis because of problems associated with detector saturation and the defini- 
tion of the time reference given by the arrival of the leading particle. The 
measurements available in the literature therefore represent the arrival time 
distributions at somewhat larger distances, where the distributions are signif- 
icantly broader. In addition, these measurements registered all charged parti- 
cles, including muons. At increasing distances from the shower core, the muons 
arrive systematically earlier than electrons and positrons and thus addition- 
ally widen the arrival time distribution. In agreement with our expectations, 
the CORSIKA-derived arrival time distribution for electrons (and positrons) 
with distances greater than 40 m from the shower axis becomes qualitatively 
similar to the parametrisation (and thus measurements). 

The differences in particle arrival time directly affect the shape of the radio 
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Fig. 7. Distribution of particle ar- 
rival times in parametrised and his- 
togrammed versions (averaged over all 
other distributions if not stated other- 
wise, normalised to a total of 10^ par- 
ticles). 



Fig. 8. Radio pulses for fully 
parametrised showers compared with 
those for showers with histogrammed 
shower evolution, energy, lateral and 
arrival time distributions. 



pulses as plotted in Fig. 8. In particular, close to the shower centre, where 
geometrical time delays are of no importance, the pulse shape directly reflects 
the particle arrival time distribution and thus becomes much narrower with 
a steeper leading edge in the histogrammed version. (Atmospheric refraction 
which is expected to slightly smear out the pulses is not taken into account 
here.) 

Additionally, the pulses amplitudes close to the shower centre are increased 
significantly with respect to those shown in Fig. 6. This is a consequence of the 
strong correlation between the energy, lateral and arrival time distributions, 
the combination of which preferably puts high-energy particles with little time 
lag close to the shower core. These particles then contribute strongly to the 
radio pulses close to the shower centre. If only one of the three distributions is 
switched back to parametrised, this amplification effect diminishes, illustrating 
how important a realistic representation of the strongly correlated particle 
distributions is. 



3.5 Particle momentum angles to the shower axis 



Finally, we additionally switch on the histogrammed particle momentum angle 
distribution. In the parametrised version, the particle momentum direction (at 
the start of a particle trajectory) was directly linked to the position of particle 
creation. The histogrammed momentum angle distribution, on the other hand, 
is strongly correlated with the particle energy distribution. 

Figure 9 shows the angular distribution of particle momenta at particle cre- 
ation relative to the shower axis. In the parametrised version, the particle 
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momenta are very collimated. The histogrammed distribution is much wider. 
Its strong correlation with particle energy is illustrated by a comparison of the 
overall distribution of electrons with that of particles with a Lorentz factor of 
60. 
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Fig. 10. Radio pulses for fully 
parametrised showers compared with 
those for fully histogrammed showers. 



Fig. 9. Particle momentum angles to 
the shower axis at particle creation (av- 
eraged over all distributions except at- 
mospheric depth if not stated other- 
wise, normalised to a total of 10^ par- 
ticles) . 



Correspondingly, and expectedly, the effect on the radio pulses is strong, as 
depicted in Fig. 10 in comparison with Fig. 8. The emission of relativistic par- 
ticles is beamed into a narrow cone along the particle momentum direction. In 
the histogrammed case, the momentum direction is much less aligned with the 
shower axis, and therefore a significant fraction of the radio signal is emitted 
at large angles to the shower axis. As a consequence, the pulses are damped 
significantly, especially in the centre region, where the high-energy particles 
with especially narrow beaming cones contribute most. 



3.6 Changes along the east-west axis 



So far we have only looked at pulses to the north of the shower centre. 
Here (and likewise to the south), switching from the parametrised to the his- 
togrammed particle distributions mainly changes the pulse amplitudes, but 
leaves the pulse shapes similar. The situation can be different to the east 
or west of the shower centre. There, the contributions from electrons and 
positrons interact in a much more complex way: the geometry under which 
the observer sees electrons and positrons, which themselves bend to the east 



and west in the geomagnetic field, becomes asymmetric. (See 12] for a detailed 
discussion.) 

In simulations with REASl, pulses at moderate to large distances in the east 
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or west from the shower centre could become bipolar, depending on the specific 
air shower configuration. The reason is that the pulses emitted by individual 
particles have a slightly bipolar structure and, if an observer starts to see only 
part of these bipolar pulses, the polarity of the overall pulse arising from the 
superposition of many individual particle pulses can change. (An example is 
shown in Fig. 18 of [12] together with a detailed discussion.) 

Such effects, however, can only occur if the air shower model (more precisely 
its particle distributions) exhibit clear "structures". This is the case for the 
parametrised distributions: the angular distribution of particle momenta is 
very narrow and the particle energy distribution has a prominent peak around 
Lorentz factors of 60. Even more pronounced bipolarities can be produced 
when the air shower model is simplified further, e.g., by neglecting the lateral 
structure of the air shower or running the simulation with long, fixed particle 
track lengths. As the CORSIKA-derived histogrammed particle distributions 
do not show any such prominent "structures" , the bipolar pulses do not occur 
in simulations of CORSIKA-histogram based showers. 
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Fig. 11. Influence of the histogrammed 
energy distribution on the radio pulses 
to the west of the shower centre. 
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Fig. 12. Influence of the histogrammed 
energy and momentum angle distribu- 
tions on the radio pulses to the west of 
the shower centre. 



Figure 11 shows the changes of pulses to the west from the shower centre 
associated with switching from the parametrised to the histogrammed energy 
distribution. The bipolarity of the 525 m pulses is washed out significantly 
because the prominent peak in the energy distribution (Fig. 3) is no longer 
present. Similarly, figure 12 demonstrates the effect of additionally switching 
on the realistic momentum angular distribution. The bipolarity of the 525 m 
pulses has diminished even more. When all of the histogrammed distribu- 
tions are switched on simultaneously, the bipolarity in the pulses vanishes 
completely. Bipolarities present in the REASl simulations can therefore be 
explained as an artifact of the simplified air shower model. 

Having illustrated in detail the effects arising from the transition to realistic, 
histogrammed particle distributions, we now investigate a second major factor 
changing the radio pulses in the transition from REASl to REAS2. 
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4 Eliminating the track length parameter 



In simulations with REASl, a free parameter, the so-called "track length pa- 
rameter" A, had to be set by the user. A denotes the average length (in g cm^^) 
of an individual particle trajectory in the simulation. It relates the longitudi- 
nal shower profile A^(^) (number of particles traversing an imaginary detector 
plane as a function of atmospheric depth X) as taken from the CORSIKA sim- 
ulations to the longitudinal particle injection profile I{X) (number of particles 
being created per dX depth interval along the shower 3jX1S clS db function of 
atmospheric depth X) via 



which ensures that the total track length integrated over all particles is con- 
stant irrespective of the value of A. A large value leads to a simulation with 
few long tracks, whereas a small value corresponds to a simulation with many 
short tracks. Even though the normalisation is correct irrespective of the value 
of A, the simulated radio signals can be influenced by the choice of a specific 
value. In particular, the angular distribution of particles evolves during their 
propagation along the trajectories. Consequently, the choice of a specific A in- 
directly influences the effective angular spread of the particle momenta. (For 
details on the role of the track-length parameter and the differences between 
fixed and exponentially distributed track lengths, see the discussion in (l2|.) 

In REASl simulations, we usually set A to a value of 36.7 g cm~^ with an 
exponential distribution of particle track lengths. As explained above, this 
implicitly sets a characteristic angular spread for the particle momenta. It is, 
however not straight-forward to motivate the choice of a particular value for 
A either from experimental data or from simulations. 

In REAS2 simulations, the situation is different. The CORSIKA-derived his- 
tograms contain detailed information on the angular spread of the particle 
momenta during the course of the full shower evolution. If A is set to inap- 
propriately large values in this case, the angular spread of particle momenta 
(and also the distribution of particle energies) starts to deviate systematically 
from the histogrammed distributions in the course of the particle trajectories. 
The simulation can then become inconsistent. 

When A is set to a small value, on the other hand, each particle always follows 
the local particle distributions well throughout its full trajectory. Trajectories 
longer than A are then effectively described by an integration over multiple, 
independent representative segments of length A. In REAS2, setting A to a 
sufficiently small value solves the problem of a particular choice for the A pa- 
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Fig. 13. Radio signal convergence in a fully histogrammed air shower when going 
from 36.7 g cm~^, exponentially distributed, to 1.0 g cm~^, fixed track lengths. 



rameter in a natural way. The radio pulses converge and become independent 
of the value of A when its value is chosen smaller and smaller (cf. Fig. 13). For 
air showers in the US standard atmosphere a value of 1 g cm~^ proves to give 
good results while the calculation is still numerically stable and efficient. 



5 Particle pairing 



In REASl, electrons and positrons are always created in pairs — a reasonable 
approximation of the situation in air showers. As explained in , part of the 



emission from electrons and positrons adds up, whereas other contributions 
cancel out. Pairwise creation of particles increases the numerical stability of 
calculating these cancellations. It is, however, well-known that air showers in 
fact exhibit a considerable negative charge excess of order ~ 20% at low elec- 
tron energy, which could influence the radio pulse strengths and polarisation 
properties. 

The histograms obtained from CORSIKA contain separate information for 
electrons and positrons, including the differences in particle multiplicity. In 
REAS2, differences in the electron and positron distributions are therefore 
taken into account in detail. To assess the importance of these differences, the 
radio emission can be calculated in REAS2 with different strategies: electrons 
and positrons can always be created in pairs ("always paired", as in REASl), 
electrons and positrons can be created in pairs (when possible), but with 
the correct ratio of electrons to positrons in each bin of phase-space ("ratio 
paired"), and electrons and positrons can be created completely independently 
( "never paired" ) with the correct ratio in each bin of phase-space. 

Figure 14 demonstrates the effect of the different "pairing modes" for ob- 
servers to the north from the shower centre. The differences are negligible. 
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Fig. 14. Differences between particle 
pairing modes in the north from the 
shower centre. 



Fig. 15. Differences between particle 
pairing modes in the west from the 
shower centre. 



For an observer to the west from the shower centre, the electron to positron 
ratio can be expected to have more influence because the geometry is asym- 
metric for the two species. A systematic difference is indeed visible between 
the "always paired" mode and the modes that correctly take into account 
the electron-to-positron ratio. Similar effects can be seen in the north-south 
polarisation component (not shown here), i.e., there is a slight effect on the po- 
larisation properties of the emission. No significant difference is visible between 
the "never paired" and "ratio paired" modes, providing additional support for 
the validity of the "small A" simulation strategy. 



6 Overall result 



So far we have only looked at isolated aspects of the air shower model and 
their influence on the calculated radio signal. Now we compare the overall 
changes between the radio emission from a fully parametrised shower (as in 
REASl, with A = 36.7 g cm~^ exponential track length distribution and 



"always paired" particles, see 12; uM) and a fully histogrammed air shower 
(REAS2, with A = 1 g cm~^ fixed track lengths and "ratio paired" particles). 

Figure 16 shows the overall difference in the calculated radio pulses for ob- 
servers along the north-south axis from the shower centre. The pulses are 
damped up to a factor of two in amplitude. Pulses close to the shower centre 
become much narrower with a much steeper rising edge. This is also reflected 
in the spectra shown in Fig. 18: close to the shower axis, the spectra are 
now much flatter. In the observing bandwidth of the LOPES experiment (40- 
80 MHz), the frequency- averaged field strength is changed only very little. 

The situation is more complex for observers along the east-west axis from the 
shower centre as shown in Fig. 17. Observers close to the shower axis receive 
a strongly damped radio pulse. For observers at larger distances the pulse 
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Fig. 16. Overall change in the simulated 
radio pulses north of the shower centre 
between REASl and REAS2. 



Fig. 17. Same as Fig. 16 but for ob- 
servers to the west of the shower centre. 



amplitude changes much less, but the bipolarity of the pulses vanishes as a 
consequence of the more realistic particle distributions. 
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Fig. 19. Same as Fig. 18 but for 
observers to the west of the shower 
centre. 



Fig. 18. Overall change in the radio sig- 
nal frequency spectra between REASl 
(thin lines) and RE AS 2 (thick lines) at 
various distances to the north of the 
shower centre. 

Figure 20 illustrates the changes in the lateral distribution of the 10 MHz 
spectral radio emission component. As expected, there is an overall damping. 
Except in the centre region, the field strength still decreases approximately ex- 
ponentially with distance, and the scale parameter of ~ 125 m for the emission 
from the histogrammed shower is not much different from that derived for the 
parametrised shower. However, a more pronounced asymmetry between north- 
south and east-west direction is observed. At 60 MHz (shown in Fig. 21) the 
overall drop in amplitude is less pronounced than at 10 MHz due to the now 
flatter frequency spectra. The scale parameter there corresponds to ~ 50 m, 
but the exponential decrease is only valid up to ~ 250-300 m distance. 

The spatial distribution of the 10 MHz radiation component and its polari- 
sation properties are also illustrated by Fig. 22. One can easily identify the 
overall damping of the emission. In addition, the radio "footprint" becomes 
more asymmetric and less circular. The polarisation characteristics remain 
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Fig. 21. Same as Fig. 20 for the 60 MHz 
frequency component. 



Fig. 20. Spectral field strength at 
10 MHz as a function of radial distance 
to the shower centre for different az- 
imuthal directions compared between 
REASl (thin lines) and REAS2 (thick 
lines). 

qualitatively very similar, even though the realistic electron to positron ratio 
is now taken into account. 
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Fig. 22. Contour plots of the parametrised (upper) and histogrammed (lower) air 
shower emission at = 10 MHz. The columns (from left to right) show the total field 
strength, the north-south polarisation component and the east-west polarisation 
component. The vertical polarisation component (not shown here) does not contain 
any significant flux. Contour levels are 0.25 /iV m~^ MHz~^ apart in E^, outermost 
contour corresponds to 0.25 //V m~^ MHz~^. White centre region has not been 
calculated. 



Figure 23 shows the changes occurring at 60 MHz, the central frequency of 
the LOPES observing bandwidth. As a consequence of the flatter spectra, 
the changes at 60 MHz are much smaller than at 10 MHz. The increasing 
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7 Pulse shape analysis 



In REAS2, detailed information on the relevant particle distributions and 
their development throughout the full air shower evolution is available. This 
information can be used to analyse how different regimes identifiable in the 
particle distributions contribute to the radio pulses. In the following we discuss 
a few examples of such analyses. 



7.1 Particle energy 



Figure 24 shows the contribution of different particle energy regimes to the 
overall radio pulse close to the shower centre. It is evident that high-energy 
particles with Lorentz factors between 100 and 1000 dominate the emission, 
although low-energy particles outnumber them by far (cf. Fig. 3). The reason 
for this is that high-energy particles radiate more efficiently than low-energy 
particles, but only into narrow beaming-cones and only from close to the 
shower axis where they are prominent. Consequently, the pulse at 525 m is 
dominated by the lower energy particles with Lorentz factors between 10 and 
100, as shown in Fig. 25. 
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Fig. 24. Contribution of different parti- 
cle energy regimes to the radio pulse at 
75 m north from the shower centre. 
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Fig. 25. Contribution of different parti- 
cle energy regimes to the radio pulse at 
525 m north from the shower centre. 



7.2 Lateral distance 



Another interesting question is how different lateral distance ranges (i.e., par- 
ticles in geometrically defined rings, independent of atmospheric depth) con- 
tribute to the emission. This is shown in Figures 26 and 27. Close to the shower 
centre, particles near the shower axis dominate the emission. This is consistent 
with the result that the high-energy particles (which are clustered close to the 
shower axis) dominate the emission close to the shower centre. Consequently, 
at larger distances, the relative importance of the particles close to the shower 
axis is lower, and lower energy particles with distances of 10-100 m from the 
axis give the dominant contribution. In all cases, particles at distances of more 
than 100 m from the axis only contribute very little to the overall pulses. 
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Fig. 26. Contribution of different lateral 
rings to the radio pulse at 75 m north 
from the shower centre. 



Fig. 27. Contribution of different lateral 
rings to the radio pulse at 525 m north 
from the shower centre. 



7.3 Atmospheric depth 



A particularly interesting question is how different stages of the air shower 
evolution contribute to the radio pulses. Figure 28 illustrates that the radio 
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emission from different slices in atmospheric depth arrive nearly simultane- 
ously at an observer close to the shower centre. The reason is that geometrical 
time-delays are small at these radial distances and that both the radio emission 
and the leading particle travel with the speed of light in the current simula- 
tions. If the refractive index profile of the atmosphere is taken into account, 
the contributions are likely to arrive at slightly different times and the overall 
pulse will get slightly broader. 

At larger distances from the shower axis, the relative arrival times between the 
emission from different air shower stages are mainly governed by geometrical 
time-delays. Figure 29 shows the radio pulse at 525 m to the north from the 
shower centre. One can clearly differentiate the contributions from the different 
air shower stages. The emission is dominated by the shower stage around the 
shower maximum (which here is at ~ 640 g cm~^) and the stage shortly 
before the shower maximum, whereas the stage after the shower maximum 
contributes less radiation. This is somewhat surprising as the particle numbers 
before and after the shower maximum are similar but the particles after the 
shower maximum are closer to the observer and thus could have been expected 
to contribute more strongly. Two effects, however, overcompensate this simple 
distance dependence: For once, the emission is strongly beamed. Emission from 
earlier shower stages is thus distributed over a larger area and therefore leads 
to stronger radio signals for observers not very close to the shower centre. 
In addition, the particles in the stage before the shower maximum traverse 
a much thinner medium than those in later stages and thus propagate for 
longer (geometrical) distances, during which they can emit radio signals. In 
general, a thinner medium leads to higher geosynchrotron radio emission per 
traversed depth in g cm~^ — if the shower is able to evolve to high particle 
multiplicities. 
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Fig. 28. Contribution of different 
shower evolution stages to the radio 
pulse at 75 m north from the shower 
centre. 
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Fig. 29. Contribution of different 
shower evolution stages to the radio 
pulse at 525 m north from the shower 
centre. 
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7.-^ Geometric height 



Figures 30 and 31 illustrate how different regions of geometrical height above 
ground contribute to the radio pulses from the air shower. The qualitative 
behaviour is similar as for the stages of atmospheric depth, but with one 
important difference: the contributions from the different height regimes con- 
tribute approximately equally to the overall pulse. This, again, is explained by 
beaming effects and the fact that radio emission scales with the geometrical 
length of the particle tracks, not the atmospheric depth traversed. 
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Fig. 30. Contribution of different 
heights to the radio pulse at 75 m north 
from the shower centre. 



50 



100 



150 



200 



250 



300 



time [ns] 

Fig. 31. Contribution of different 
heights to the radio pulse at 525 m 
north from the shower centre. 



8 Conclusions 



We have analysed in detail how improving the air shower model from purely 
parametrised particle distributions to realistic histograms obtained with COR- 
SIKA changes the radio emission predicted for a vertical 10^^ eV proton- 
induced reference shower. The approach of keeping the well-understood radio 
code and only substituting the air shower model is ideally suited to understand 
the changes in detail. 

The effects of the transition to realistic, histogrammed particle distributions 
have been evaluated step by step. One important change is introduced by the 
transition to the histogrammed arrival time distribution, which produces sig- 
nificantly narrower pulses and, consequently, flatter frequency spectra, in the 
shower centre region. (Measurements of radio pulse shapes close to the shower 
core in comparison with the simulations presented here could therefore be used 
to study the particle arrival time distributions in air showers.) As expected, by 
far the most important difference is caused by taking into account the realistic, 
much broader distribution of particle momentum angles to the shower axis. 
The histogrammed distribution dampens the radio emission significantly, as 
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it causes a considerable fraction of the radiation to be emitted at large angles 
to the shower axis. It also became clear that the strong correlations between 
the different (e.g., energy, lateral and arrival time) distributions play an im- 
portant role in shaping the radio emission and cannot be neglected. Another 
major improvement is the fact that the adoption of longitudinal shower evolu- 
tion profiles obtained from CORSIKA allows us to analyse shower to shower 
fluctuations and differences related to the mass of the primary particle. 

As a consequence of the much more realistic distributions in general, a qual- 
itative change in the pulse shapes measured by observers along the east-west 
axis from the shower centre arises. Bipolarities in the pulses that were present 
in the simulations with parametrised air showers vanish when the realistic 
distributions are adopted. The bipolarities present in our earlier simulations 
seem to be artifacts of the over-simplified air shower model. 

The track length parameter A was eliminated in a natural way by simulating 
long tracks with an appropriate number of independent, representative short 
tracks. This approach ensures that the particle distributions throughout the 
tracks always represent the local distributions and that ionisation losses along 
a track can be safely neglected. The radio pulses are shghtly enhanced by 
this choice of A, compensating for some of the damping introduced by the 
histogrammed air shower model. 

Different modes of creating particles independently or in pairs, with equal 
numbers of electrons and positrons or the correct ratio, have been investigated. 
Taking into account the correct electron to positron ratio only slightly changes 
the emission along the east-west axis from the shower centre. 

Another change is visible in the azimuthal distribution of the radio emission. 
The radio "footprint" becomes significantly more asymmetric, as the damping 
of the radio signal is more pronounced in the east and west than in the north 
and south. The radial decrease of the radio emission, however, stays approx- 
imately exponential, and the scale parameter does not change significantly 
either for the vertical shower studied here. 

The new code allows us to analyse how different regimes identifiable in the 
shower particle distributions contribute to the radio pulses measured on the 
ground. In particular, different stages of the air shower evolution contribute to 
different time-windows of the radio pulses at moderate to large distances from 
the shower centre, where geometric time-delays play an important role. This 
shows that, in principle, information on the air shower evolution (such as the 
position of the shower maximum, which in turn is related to the mass of the 
primary particle) is encoded in the radio signal and could be accessible through 
radio measurements of cosmic ray air showers. The contributions of different 
slices in geometrical height above ground illustrate that beaming effects play 
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an important role and that the geosynchrotron radio emission scales with the 
geometrical length of particle tracks, not the atmospheric depth traversed. 
This means that the same number of particles in a less dense medium produces 
more geosynchrotron radiation than in a dense medium. The analyses of the 
energy and radial distance regimes demonstrate that in the shower centre, 
high-energy particles located close to the shower axis dominate the emission, 
whereas at distances of a few hundred metres, lower-energy particles located 
further away from the shower axis dominate. 

With REAS2, a next-generation, well-tested Monte Carlo code for the calcula- 
tion of geosynchrotron radio emission from cosmic ray air showers is available. 
As a next step, we will repeat our earlier analysis of the influence of specific air 
shower parameters on the associated radio emission (cf. [l^) with radio simu- 
lations based on CORSIKA-simulated showers. Additionally, we will assess the 
role of shower-to-shower fluctuations and the primary particle composition. In 
the future, we will incorporate the effect of the atmospheric refractive index 
proflle in combination with the presence of systematic lateral shifts between 
the electron and positron distributions to take into account Cerenkov-like ra- 
diation contributions. 



Acknowledgements 



The authors would like to thank T. Pierog and D. Heck for their outstand- 
ing support regarding many CORSIKA-related aspects of this work and the 
anonymous referee for several valuable comments that helped to improve the 
article. T.H. acknowledges very helpful discussions with S. Buitink, A. Haungs 
and H. Falcke. R.E. would like to thank N. Kalmykov and A. Konstantinov 
for many useful discussions. 



References 

[1] H. Falcke, W. D. Apel, A. F. Badea, et al.. Detection and imaging of 
atmospheric radio flashes from cosmic ray air showers. Nature 435 (2005) 
313-316. 

[2] A. Horneffer, W. D. Apel, A. F. Badea, et al.. Radio Detection of Cos- 
mic Rays with LOPES, in: Int. J.Mod.Phys. A; Proc. of the ARENA2005 
conference, DESY Zeuthen, Germany, no. 21S1, World Scientiflc, 2006, 
pp. 168-181. 

[3] W. D. Apel, T. Asch, A. F. Badea, et al.. Progress in air shower radio 
measurements: Detection of distant events, Astropart. Physics 26 (2006) 
332-340, ,astro-ph/0607495. 



23 



[4] T. Huege, W. Apel, T. Asch, A. Badea, et al., Radio detection of cos- 
mic ray air showers with LOPES, in: Nuclear Physics B, Proceedings 
Supplements; Proc. of the CRIS2006 conference, Catania, Italy, 2006, 
'astro-p h/0609427 . 

[5] D. Ardouin, A. Belletoile, D. Charrier, et al.. Radio-detection signature 
of high-energy cosmic rays by the CODALEMA experiment, Nucl. Instr. 
Meth. A 555 (2005) 148-163. 

[6] D. Ardouin, A. Belletoile, D. Charrier, R. Dallier, L. Denis, P. Eschstruth, 
T. Gousset, et al., Radioelectric field features of extensive air showers 
observed with CODALEMA, Astropart. Physics 26 (2006) 341-350. 

[7] H. R. Allan, Radio Emission From Extensive Air Showers, Prog, in Ele- 
ment, part, and Cos. Ray Phys. Vol. 10 (1971) 171-302. 

[8] H. Falcke, P. W. Gorham, Detecting radio emission from cosmic ray air 
showers and neutrinos with a digital radio telescope, Astropart. Physics 
19 (2003) 477-494. 

[9] D. A. Suprun, P. W. Gorham, J. L. Rosner, Synchrotron radiation at radio 
frequencies from cosmic ray air showers, Astropart. Physics 20 (2003) 
157-168. 

[10] T. Huege, H. Falcke, Radio emission from cosmic ray air showers. Co- 
herent geosynchrotron radiation. Astronomy & Astrophysics 412 (2003) 
19-34. 

[11] Q. Luo, Coherent synchrotron emission from cosmic ray air showers. 
Monthly Not. Roy. Astronom. Soc. 370 (2006) 2071-2078. 

[12] T. Huege, H. Falcke, Radio emission from cosmic ray air showers. Monte 
Carlo simulations. Astronomy & Astrophysics 430 (2005) 779-798. 

[13] T. Huege, H. Falcke, Radio emission from cosmic ray air showers: Simu- 
lation results and parametrization, Astropart. Phys. 24 (2005) 116. 

[14] S. J. Sciutto, Aires: A system for air shower simulations (version 2.2.0), 
astro-ph/9911331, see also http://www.fisica.unlp.edu.ar/auger/aires, 

[15] M. A. Duvernois, et al., Geosynchrotron radio pulse emission from exten- 
sive air showers: Simulations with AIRES, in: Proc. of the 29th ICRC, 
Pune, India, 2005, pp. 311-+. 

[16] http:/ /www. irs.inms.nrc.ca/EGSnrc/EGSnrc. html. 

[17] R. Engel, N. N. Kalmykov, A. A. Konstantinov, Simulation of Radio 
Signals from 1-10 TeV Air Showers using EGSnrc, in: Int. J. Mod. Phys. 
A; Proc. of the ARENA2005 conference, DESY Zeuthen, Germany, no. 
21S1, World Scientific, 2006, pp. 65-69. 

[18] D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, T. Thouw, CORSIKA: 
A Monte Carlo Code to Simulate Extensive Air Showers, FZKA Report 
6019, Forschungszentrum Karlsruhe (1998). 

[19] F. Nerling, J. Bliimer, R. Engel, M. Risse, Universality of electron distri- 
butions in high-energy air showers — Description of Cherenkov light pro- 
duction, Astropart. Physics 24 (2006) 421-437. 

[20] M. Kobal, Pierre Auger Collaboration, A thinning method using weight 
limitation for air-shower simulations, Astroparticle Physics 15 (2001) 259- 



24 



273. 

[21] S. Ostapchenko, Non-linear screening effects in high energy hadronic in- 
teractions, Phys. Rev. D74 (2006) 014026. 

[22] S. Ostapchenko, On the re-summation of enhanced pomeron diagrams, 
Phys. Lett. B636 (2006) 40-45. 

[23] S. A. Bass, et al. Prog. Part. NucL Phys. 41 (1998) 225. 

[24] M. Bleicher, et al., J. Phys. G 25 (1999) 1859. 

[25] W. R. Nelson, H. Hirayama, D. W. O. Rogers, The EGS4 Code 
System, SLAG Report 265, Stanford Linear Accelerator Genter, 
http:// www.slac.stanford.edu/pubs/slacreports/slac-r-265.html (1985). 

[26] K. Greisen, Prog. Gosmic Ray Phys. Vol. 3 (1956) 1. 

[27] K. Kamata, J. Nishimura, Prog. Theoret. Phys. Suppl. Vol. 6 (1958) 93. 

[28] G. Agnetta, M. Ambrosio, G. Aramo, et al.. Time structure of the exten- 
sive air shower front, Astropart. Physics 6 (1997) 301-312. 



25 



